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Abstract 

In this article we address the regularization of the ill-posed problem of determining 
the local volatility surface (as a function of time to maturity and price) from market 
given option prices. We integrate the ever-increasing flow of option price information 
into the well-accepted local volatility model of Dupire. This leads to considering both 
the local volatility surfaces and their corresponding prices as indexed by the observed 
underlying stock price as time goes by in appropriate function spaces. 

The parameter to data map consists of a nonlinear operator that maps the (vari- 
able) diffusion coefficient of a parabolic initial value problem into its solutions evalu- 
ated at certain sets. We tackle the inverse problem by convex regularization techniques 
in appropriate Bochner-Sobolev spaces. 

As a preparation, we prove key regularity properties that enable us to apply convex 
regularization techniques. This forward framework is then used to build a calibration 
technique that combines online methods with convex Tikhonov regularization tools. 
Such procedure is used to solve the inverse problem of local volatility identification. As 
a result, we prove convergence rates with respect to noise and a corresponding Morozov 
discrepancy principle for the regularization parameter. We conclude by illustrating and 
validating the theoretical results by means of numerical tests with synthetic as well as 
real data. 

Keywords: Local Volatility Calibration, Convex Regularization, Online Estimation, Mo- 
rozov's Principle, Convergence Rates. 

Introduction 

A number of interesting problems in nonlinear analysis are motivated by questions from 
mathematical finance. Among those problems, the robust identification of the variable 
diffusion coefficient that appears in Dupire's local volatility model [7, 12] presents substan- 
tial difficulties for its nonlinearity and ill-posedness. In previous works tools from Convex 
Analysis and Inverse Problem theory have been used to address this problem. See [6] and 
references therein. 
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In this work, we incorporate the fact that as time evolves more data is available for 
the identification of Dupire's volatility surface. Thus we develop an online approach to the 
ill-posed problem of the local volatility surface calibration. Such surface is characterized by 
a non-negative two- variable function a = a(r, K) of the time to expiration r and the strike 
price K. 

In what follows, we consider that the local volatility surfaces are indexed by the observed 
underlying asset price So- The reason for that stems from the fact that if we try to use 
information of prices observed on different dates, there is no financial or economical reason 
for the volatility surface to stay exactly the same. Thus, in principle we may have different 
volatility surfaces, although such change may be small. 

Let us quickly review the standard Black-Scholes setting and Dupire's local volatility 
model. Recall that an option or derivative is a contract whose value depends on the value 
of an underlying stock or index. Perhaps the most well known derivative is an European 
call option, where the holder has the right (but not the obligation) to buy the underlying 
at time t = T for a strike value K. We shall denote the stochastic process defining such 
underlying S(t) = S(t, u), where as usual we assume that it is an adapted stochastic process 
on a suitable filtered probability space (Q, % ,W, P), where F = {Fi} t£ R is a filtration [15]. 

It is well known [7, 12, 15] that the value C of an European call option with strike K 
and expiration T = t + r, where t is the current time, satisfies: 

C(r = 0,K) = (S -K) + , forK>0, () 
lim C(t, K) = 0, for r > 0, U 

K— >+oo 

lim C(t,K) = S , for r>0 

where b is the difference between the continuously compounded interest and dividend rates 
of the underlying asset. In what follows, we assume that such quantities are constant. 
Defining the diffusion parameter a(r,K) = <j(t, K) 2 /2, Problem (1) leads to the following 
parameter to solution map: 

F : D(F) CX — > Y 

a G D(F) i— )■ F(a) = C eY 

where X and Y are Hilbert spaces to be properly defined below. D(F) is the domain of the 
parameter to solution map (not necessarily dense in X) and C = C(a, r, K) is the solution 
of Problem (1) with diffusion parameter a. 

The inverse problem of local volatility calibration, as it was tackled in previous works 
[4, 5, 6, 8], consists in given option prices C , find an element a of D(F) such that F(a) = C . 
The operator F is compact and weakly closed. Thus, this inverse problem is ill-posed. In 
[4, 5, 6, 8] different aspects of Tikhonov regularization were analyzed. It is characterized by 
the following: Find an element of 

argmin { ||F(a) — C||y + a/ ao (a)} subject to a E D(F) C X, 

where f ao is a weak lower semi-continuous convex coercive functional. The analysis presented 
in [4, 5, 6, 8] was based on an a priori choice of the regularization parameter with convex 
regularization tools. 

In contrast, in the present work we explore the dependence of the local volatility surface 
on the observed asset price in order to incorporate different option price surfaces in the same 
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procedure of Tikhonov regularization. More precisely, we consider the map 



u : d(u) c x i — ► y 

A e D{U) .— ► U{A):S e[S min ,S max ]^C(S,a(S)) 

where C(S, A(S)) is the solution of (1) with S = S and a 2 /2 = a(S). Moreover, A maps 

5 G [S min , S m3X \ to a(S) G D(F) in a well-behaved way. 

In this context the inverse problem becomes the following: Given a family of option 
prices C G y, find A G D(U) such that U(A) = C. We shall see that the operator U is 
also compact and weakly closed. Thus, this problem is also ill-posed. The corresponding 
regularized problem is defined by the following: 

Find an element of 

/ \\F(a(S)) -C(S) \\ydS + af M (A) \ subject to A G D(U). 
The main contributions of the current work are the following: 

Firstly, we extend the local volatility calibration problem to local volatility families. 
This new setting allows incorporating more data into the calibration problem, leading to an 
online Tikhonov regularization. We prove that the so-called direct problem is well-posed, 
i.e., the forward operator satisfies key regularity properties. This framework generalizes 
in a nontrivial way the structure used in previous works [4, 5, 6, 8] since it requires the 
introduction of more tools, in particular that of Bochner spaces. 

Secondly, in this setting, we develop a convergence analysis in a general context, based 
on convex regularization tools. See [19]. 

Thirdly, we establish a relaxed version of Morozov's discrepancy principle with conver- 
gence rates. This allows us to find the regularization parameter appropriately for the present 
problem. See [2, 17]. 

The paper is divided as follows: 

In Section 1, we present the setting of the direct problem. In Section 2, we define 
properly the forward operator and prove some key regularity properties that are important 
in the analysis of the inverse problem. This is done in Theorem 1 and Propositions 4, 5, 

6 and 7. In Section 3, we tie up the inverse problem with convex Tikhonov regularization 
under an a priori choice of the regularization parameter. This is done in Theorems 2, 3, 4 
and 5. In Section 4 we establish the Morozov discrepancy principle for the present problem 
with convergence rates. This is done in Theorems 7 and 8. Illustrative numerical tests are 
presented in Section 5. 



1 Preliminaries 

We start by setting the so-called direct problem. It is based on the pricing of European call 
options by a generalization of Black-Scholes-Merton model. 

Performing the change of variables y := \og(K / So) and r := T — t on Problem (1) and 
defining 

u(S , t, y) := C(S , r + t, S e y ) and a(S , r, y) := ^a 2 (S , r + t, S e y ), 
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it follows that u(So,r,y) satisfies 



d 2 u 
dy 2 



dy) dy 
u(t = 0,y) 
lim u(r,y) 

y— >+oo 

lim u(r,y) 

y— >— oo 



t > 0, y e 

S {l-e y ) + , for yeR, 

0, for r > 0, 

Sq, for r > 0. 



(2) 



Note that, a and a are assumed strictly positive and are related by a smooth bijection 
(since a > 0). Thus, in what follows we shall work only with the local variance a instead of 
volatility a. This simplifies the direct and inverse problems analysis. 

Let I C R be a possibly unbounded open interval, denote by D = (0, T)x I the set where 
problem (2) is defined. From [8] we know that (2) has a unique solution in W 2 f oc (D), the 
space of functions u : (r, y) G D \-t u{r,y) G E such that, it has locally squared integrable 
weak derivatives up to order one in r and up to order two in y. 

For latter purposes, we define the set where the diffusion parameter a lives. For fixed 
e > 0, take scalar constants a\,d2 G K such that < a± < ci2 < +oo and a fixed function 
do G H 1+£ (D), with ao < a < di- Define 



g := {a G a + H 1+£ (D) : a x < a < a 2 } 



(3) 



Note that Q is weakly closed and has nonempty interior under the standard topology of 
H 1+£ (D). See the first two chapters of [5, 6] and references therein. 



2 The Forward Operator 

Since we assume that the local variance surface is dependent on the current price, we have 
to introduce proper spaces for the analysis of the problem. As it turns out, we have to 
make use of Bochner spaces techniques. For general results on the Bochner integral see 
[11, 18, 23]. The main reference for this section is [14]. 

We start with some definitions. Given a time interval, say [0,T], the realized prices 
S(t) vary within [S m m : S m «x\- For technical reasons we perform the change of variables 
s = S(t) — S min and denote S = S max — S min . Thus s G [0, S). Hence, for each s, we denote 
a(s) := a(s,r,y) the local variance surface correspondent to s. 

Definition 1. Given A G L 2 (0 , S , H 1+e (D)) , with i : s 4 a(s) (see [23]), we define its 
Fourier series A = {d{k)}k^i by 

I rS 1 rO 

a(k) := — / a(s) exp(— iksir/S)ds + — / a(— s) exp(— iksix/S)ds. 
2S J 2b J_ s 

It is well defined, since {s y a(s) exp(— ihs2ir/S)} is weakly measurable and 

L 2 (0,S,H 1+£ (D)) C L^O.S'.ff^p)) 

by the Cauchy-Schwartz inequality. 

Now we define a class of Bochner-type Sobolev spaces: 
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Definition 2. Let H l {0, S, H 1+£ (D)) be the space of A E L 2 (0, S, H 1+£ {D)) , such that 

\\A\l, := £(1 + \k\ l na(k)f Hl+HD)c <oc 1 

k&L 

where H 1+£ (D)c = H 1+£ (D) © iH 1+£ (D) is the complexification of H 1+£ (D). Moreover, 
H l (0, S, H 1+£ (D)) is a Hilbert space with the inner product 

(A, A) i := ^(1 + \k\ l ) 2 {a(k),~a(k)) m+ . (D)c . 

Proposition 1. [14, Lemma 3.2] Fori > 1/2, each A E H l (0, S, H 1+£ (D)) has a continuous 
representative and the mapping i\ : H l (0, S, H 1+£ (D)) C(Q, S, H 1+£ (D)) is continuous 
(bounded). Moreover, we have the estimate 

/ oo 1 \ 3 
s^Ju( S )\\ mt , m < Ml ^^ Tw j . (4) 

Defining (A, x) H x+e {D) := {s ^ (a(s),x)} for x E H 1+£ (D) and A E H l (0, S, H l+£ (D)) , we 
have that (A,x) H i+e^ E H l ([0,S]) and \\(A,x) h ^+^(d)\\h'([o,s]) — \\A\\i\\x\\ H i+e^. More- 
over, for every A,B E L 2 (0, S, H 1+£ (D)), 

{A,B) L i(p, s ,m+z{D)) = ^2(a(k),b(k)) H i+. {D)c . 

fcez 



In what follows, we shall need the following technical lemma: 

Lemma 1. Assume that I > 1/2. If the sequence {A n } n en converges weakly to A in 
H l (0, S, H 1+£ (D)), then, the sequence {a k {s)}km weakly converge to a(s) in H 1+£ (D) for 
every s E [0, S]. 

Proof: Take a {^4 n } ng N and A as above. We want to show that, given a weak zero neigh- 
borhood U of H 1+£ (D), then for a sufficiently large n, a n (s) — a(s) E U for every s E [0, S}. 
A weak zero neighborhood U of H 1+£ (D) is defined by a set of a%, ...,ax E H 1+£ (D) and 
an e > such that g E H 1+£ (D) is an element of U if max^^...^ \ (g, a n )\ < e. 

Since the immersion H\[0,S]) C([0, B\) is compact and H\[0,S]) is reflexive, it 
follows that each weak zero neighborhood of H l ([0, S]) is a zero neighborhood of C([0, S]). 
Furthermore, from Proposition 1 we know that (A,at)H 1 + e (D) ^ H l ([0,S]) with its norm 
bounded by ||.A||z||a||#i+e(£>), for every n E N and a E H 1+£ (D). Thus, we take the smallest 

closed ball centered at zero, B, which contains {A,a k )m+ s {D) with k = 1,...,K and every 
(An, o;fc)ij'i+e(x)) with n E N and k = 1,...,K. Therefore, choosing e > as above, it is 
true that for each k = 1,...,K, there are fk,i, fk,M(k) e H l ([Q,S]) and r\ k > 0, such 
that ||/||c([o,5]) < e for every f E B with max m= i i ... jM(fc) \{f,f k ,m}\ < Vk- Hence, we define 
Ck, m ■= ock <8> fk, m E H l (0, S, H 1+£ (D))* and the weak zero neighborhood A = n k=1 A k of 
H l (Q,3,H 1+e (D)) with 

A k :={AEH l (0,S,H 1+£ (D)) : \(A,C k , m )\ < Vk , m = 1, M(k)}. 

As A is a weak zero neighborhood of H l (0, S, H l+£ (D)), it is true that for sufficiently large 
n, A n — A E A, which implies that a n (s) — a(s) E U for every s E [0,5], i.e., {a n (s)} n( =N 
weakly converges to a(s) for every s E [0, S). □ 



5 



Define: 

£2 := {A G H l (0, S, H l+£ (D)) : a(s) G Q, Vs G [0, S}}, (5) 

i.e., each .4. in £3 is the map A : s G [0, S] (-> a(s) G Q. Note that £3 is the space of Q-valued 
paths, with Q defined in (3). 

Proposition 2. For I > 1/2, the set £3 is weakly closed and its interior is nonempty in 
H l (0,S,H 1+£ (D)). 

Proof: By Lemma 1 and the fact that Q is weakly closed it follows that £3 is weakly closed. 
The assertion that the interior of £3 is nonempty follows from the fact that 

H l (0, S, H l+£ (D)) C(0, S, H 1+£ (D)) 

is continuous and bounded (note that, given e > 0, A — {s (->■ a(s)} with a + e < a(s) < a + e 
for every s G [0, S] is in the interior of £3). □ 

We stress that, in what follows, we always assume that / > 1/2, since it is enough to 
state our results concerning regularity aspects of the forward operator. 

We define below the forward operator, that associates each family of local variance 
surfaces to the correspondent family of option price surfaces, determined by (2). Thus, for 
a given a G Q we define: 

W:£3 — ► L 2 (0,S,W*' 2 {D)), 

A i— >• U{A) : s G [0, S\ h> F(s, a(s)) G W*' 2 (D), 

where [U(A)](s) = F(s,a(s)) := u(s,a(s)) — u(s,a ) and u(s,a) is the solution of (2) 
with local variance a. The following results state some regularity properties concerning the 
forward operator. 

Proposition 3. The operator F : [0, S] X Q — > W 2 ,2 (D) is continuous and compact. 
Moreover, it is sequentially weakly continuous and weakly closed. 

We define below the concept of Frechet equi-differentiability for a family of operators. 

Definition 3. We call a family of operators {J-" s : Q — > Wl' 2 (D) \ s G [0,>S]} Frechet 
equi-differentiable, if for all a G Q and e > 0, there is a 5 > 0, such that 

sup \\T t (a + h)- T s (a) - F' s (a)h\\ < e\\h\\, 

s£[Q,5] 

for ||/i||^i+e(d) < $ an d ^(a) the Frechet derivative of J~ s (-) at a. 
Thus, we have the following proposition. 

Proposition 4. The family of operators {F(s, ■) : Q — > Wl' 2 (D) \ s G [0, S] } is Frechet 
equi-differentiable. 

Proof: Given oGQ and e > 0, define w = F(s, a + h) — F(s, a) — d a F(s, a)h, it is equivalent 
to w — u(s, a + h) — u(s, a) — d a u(s, a)h. We denote v := u(s, a + h) — u(s, a). Thus, by 
linearity w satisfies 

— W T + ^(Wyy - Wy) + OW y = fl(V ' yy ~ Vy), 

with homogeneous boundary condition. Such problem does not depend on s, as a is inde- 
pendent of s. From the proof of Proposition 3 (see also [8]), we have 

|M| W 1,2 (D) < C\\h\\ L 2 {D) \\v\\ w l,2 {D) 
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By the continuity of the operator F, given e > we can chose h G H 1+£ (D) with 



, such that IMIh/^^d) < e /C an d thus the assertion follows. □ 
The following theorem is the principal result of this section, since it states some properties 
that are at the core of the inverse problems analysis [10, 19]. For its proof see Appendix 
A.3. 

Theorem 1. The forward operator U : H — > L 2 (0, 5, Wl ,2 (D)) is well defined, continuous 
and compact. Moreover, it is sequentially weakly continuous and weakly closed. 

The next results state necessary conditions for the convergence analysis. See [10, 19]. 
Its proof is in the Appendix A. 3. 

Proposition 5. The operator U(-) admits a one sided derivative at A £ Q. in the direction 
H, such that A + H G H. The derivative U'(A) satisfies 



U\A)H 

Moreover, U'(A) satisfies the Lipschitz condition 



L 2 (0,S,W^' 2 (D)) v v " 



U'{A)-U'{A + U) 



for all A, H G H such that A, A + H G H. 

The following result is a consequence of the compactness of W(-). 
Proposition 6. The Frechet derivative of the operator U(-) is injective and compact. 

Proof: Take H G ker (u'(A)^. Thus, from the proof of Proposition 5, we have 

h(s) ■ (Uyy - Uy) = 0. 

However, for each t } G = the solution of 

9 T G = \ [ply - d y ) (a(s)G + bG) 
G | r=0 = 8{y) , 

i.e., G is the Green's function of the Cauchy problem above. Thus, G > for every y,r > 
and s G [0,5*]. Therefore h{t) = 0. Since this holds for every s G [0,5], then the result 
follows. □ 
We now make use of the bounded embedding of L 2 (0, 5, W^{D)) into L 2 (0, 5, L 2 (D)), 
since it implies that U satisfies the same results presented above with L 2 (0, 5, L 2 (D)) instead 
of L 2 (0, 5, Wl' 2 {D)). Thus, we characterize the range oiU'{A) as a subset of L 2 (0, 5, L 2 {D)) 
and the range of U'(A)* as a subset of H l (0, S, H 1+£ (D)) in order to proceed in Section 3 
the convergence analysis. 

Proposition 7. The operator W (A^)* has a trivial kernel. 
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Proof: For simplicity take b = 0. Denote 

£ := -d T + a{dyy - d y ) 

the parabolic operator of Equation (2) with homogeneous boundary condition and Q Uyy -u y 
the multiplication operator by u yy — u y . Thus, for each s G [0,S], we have d a u(s,a(s)) = 
£~ 1 G U -u. , where C~ l is the left inverse of £, with null boundary conditions. By definition 
of 

U'{A)* : L 2 (0, S, L 2 {D)) — > H l (0, S, H 1+£ (D)), 

we have, 



VHG H l (0,S,H 1+£ (D)) and V Z G L 2 (0,S,L\D)), with $ = W(.A)*2. Thus, given any 
Z G ker (u'(A)*\ it follows that 

= (U'(A)H, Z) LH0AL2{D)) = [ {£~ l Gu yy - Uy h{s), z{s)) L2{D) ds 



{Gu yy ~u y h(s), [C ^ z(s)) L2(D) ds = (g Uyy „ Uy h(s),g(s)) L2{D) ds 

where g is a solution of the adjoint equation 

9t + (ag) yy + (ag) y = z 

for each s G [0,5*], with homogeneous boundary conditions. Since z(t) G L 2 (D), we have 
that g(s) G H 1+£ (D) (see [16]) and g G L 2 (0, S, H 1+£ (D)). Since Q > 0, from the proof of 
Proposition 6 and the fact that h G H l (0, S, H 1+£ (D)) is arbitrary, it follows that g = 0. 
Therefore Z = almost everywhere in s G [0, 5]. It yields that ker (U'(a)*) = {0}. □ 

Remark 1. From the last proposition it follows that 



ker{W'(^)} = {0} =>n{(u'(A)y} = H l (0,S,H 1+£ (D)). 

In other words, the range of the adjoint operator of the Frechet derivative of the forward 
operator U at A is dense in H l (0, S, H 1+£ (D)) . 

To finish this section we shall present below the tangential cone condition for U. It 
follows almost directly by the above results and Theorem 1.4.2 from [5]. 

Proposition 8. The map U(-) satisfies the local tangential cone condition 

U(A) -U(A) -U'(A)(A- A) <7 U(A)-U(A) (6) 

for all A, A in a ball B(A*,p) C H with some p > and 7 < 1/2. 

As a corollary we have the following result: 
Corollary 1. The operator U is injective. 
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3 The Inverse Problem 



Following the notation of Section 2, we want to define a precise and robust way of relating 
each family of European option price surfaces to the corresponding family of local volatility 
surfaces, both parameterized by the underlying stock price. We first present an analy- 
sis of existence and stability of regularized solutions, then we establish some convergence 
rates. We also prove Morozov's discrepancy principle for the present problem with the same 
convergence rates. 

3.1 The Regularized Problem 

The inverse problem of local volatility calibration can be restated as: 

Given a family of of European call option price surfaces U = {s H- u(s)} G L 2 (0, S, L 2 (D)), 

find the correspondent family of local variance surfaces i = {s4 G H, satisfying 

U = U(A). (7) 

We call U the observable variable and A the unknown. This is an idealized situation since 
the model (7) has no uncertainties associated. Thus, to be more realistic, we assume that 
we can only observe corrupted data U s , satisfying a perturbed version of (7), 

U & =U + £ = U{A) + £ (8) 

where £ = {t i— > E(t)} compiles all the uncertainties associated to this problem and U is the 
unobservable noiseless data. We assume further that, the norm of £ is bounded by the noise 
level 5 > 0. Moreover, for each t G [0, T], we assume that < 5/T. These hypotheses 

imply that 



\\U 5 -U\\ L 2 [oWD)) < 5 and \\u\t) - u(t) \\ L 2 (D) < 5/T for every t G [0,T]. (9) 

Proposition 1 gives that U(-) is compact, implying that the associated inverse problem 
is ill-posed. It means that such inverse problem cannot be solved directly in a stable way. 
Hence, we must apply regularization techniques. This, roughly speaking, relies on stating 
the original problem under a more robust setting. More specifically, instead of looking for 
an A s G H satisfying (8), we shall search for an A s G H minimizing the following Tikhonov 
functional 

?%M) = \\U S -U(A)\\l 2mL2m + af Ao (A)- (10) 

The functional f^ Q , which will be made more specific later, incorporates "a priori" 
information that we may have on Ao and has the goal of stabilizing the problem. 

We shall see later that, under this setting, such minimizers are, in an appropriate sense, 
good approximations for the solution of (7). We recall that the forward operator is injective. 

3.2 Some Properties of Minimizers 

In order to guarantee the existence of stable minimizers for the functional (10), we assume 
that /_4 : — > [0, oo] is convex, coercive and weakly lower semi-continuous. A classical 
reference on convex analysis is [9]. Note that, these assumptions are not too restrictive, 
since they are fulfilled by a large class of functionals on H l (0, S, H 1+e (D)). A canonical 



9 



example is fj, (A) = \\A — -4q||#!( sh i+ s r D \\, which is leads us to the classical Tikhonov 
regularization. 

Recall that U is weakly continuous and H is weakly closed. Combining that with the 
required properties of f Ao , we have that 

^ a (M) = {Ae£l\T 5 A{ha (A)<M} 

are weakly pre-compact and the restriction of U to /z Q (M) is weakly continuous. 
As a consequence of this, we have the following three theorems. See [19]. 

Theorem 2 (Existence). For every U 5 G L 2 (0, S, L 2 (D)) , there exists at least one element 
of H minimizing a (-), the functional defined in (10). 

Before stating the next theorem, we need the following definition: 

Definition 4 (Stability). If A is a minimizer of (10) with data U, then it is called stable 
if for every sequence {£4}fceN C L 2 (0, S, W 2 ' 2 (D)) converging strongly to U, the sequence 
{Ak}ken C H of minimizers of J*j^ a (-) has a subsequence converging weakly to A. 

Theorem 3 (Stability). The minimizers of (10) are stable in the sense of Definition 4- 

The following result states that, when the noise level 5 and the regularization parameter 
a = a(5) vanish, then we can find a sequence of minimizers of (10) converging weakly to 
the solution of (7). 

Theorem 4 (Convergence). Assume that (7) has a solution and the map a : (0, oo) — > 

(0, oo), satisfies 

5 2 

lim a (5) = and lim — — = 0. (11) 

<5^o s-^o a(6) 

Moreover, we assume further that the sequence {Sk}ken converges to and the elements of 
{Uk}kem, with Uk = U 5k , satisfy \\U — Uk\\ < 5k, with U the noiseless data of (7). Denote 
a(Sk) by ak for every k G N. Then, every sequence {Ak}ken of minimizers of J-^ (■), 
converges weakly to A\ the unique solution of (7), with fj^{Ak) —> /a (^)- 

Therefore we can conclude that the present regularization procedure gives us reliable 
approximations for the solutions of Problem (7). Note that, such approximations depend 
on the magnitude of 5 and thus on the choice of a and f Ao . Thus, Theorem 4 says the 
smaller 5 is, if a is properly chosen, the less dependent on regularization the solutions are. 

3.3 A Convergence Analysis 

Making use of convex regularization tools, we provide some convergence rates with respect 
to the noise level. Thus, we need some abstract concepts, as the Bregman distance related 
to /U , g-coerciveness and the source condition related to operator U. Such ideas were 
also used in [4, 5, 6, 8], but here they are extended to the context of online local volatility 
calibration. For the definitions of Bregman Distance and g-coerciveness see Appendix A.l. 

In what follows we always assume that (7) has a (unique) solution which is an element 
of the Bregman domain T>B{fA )- 

Before stating the result about convergence rates, we need the following two auxiliary 
lemmas. The first one introduces the so-called source condition. See [19]. 
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Lemma 2. For every £t G dfj^i^AJ), there exists G L 2 (0, S, L 2 (D)) and 

£ G H\0,S, H l+£ (D)) 

such that 

e=[U\A)XJ + <g (12) 
holds. Moreover, can be taken such that H^Hi^o.s.Hi+^D)) arbitrarily small. 

Lemma 2 follows from the fact that TZ{U'{A ] )*) is dense in H l (0, S, H l+£ (D)). See Propo- 
sition 7 in Section 2. Observe also that, we identify L 2 (0, S, L 2 (D))* and H l (0, S, H 1+£ (D))* 
with L 2 (0, S, L 2 (D)) and H l (0, S, H 1+£ (D)), respectively, since they are Hilbert spaces. 

Now, the main result of this section. 

Theorem 5 (Convergence Rates). Under the hypotheses of this section, let the map a : 
(0,oo) — > (0, oo) be such that a(5) ~ 5. Furthermore, assume that ,/U (-) q-coercive with 
constant (, with respect to the norm of H l (0, S, H 1+£ (D)). Then 

D e (Al,A ] ) = 0(8) and \\U(A s a ) -U s \\ = 0(6). 



Proof: Let A^ and A s a denote the solution of (7) and the minimizer of (10). It follows that, 

\\U(A s a ) - U s \\ 2 + af M {A s a ) < \\U{A) -U 6 \\ 2 + af Ao (A ] ) < $ 2 + af Ao (A ] )- 
Since, D v (Ai,A*) = f Ao {K) ~ f^i) ~ " ^ h J Lemma 2 we have that ' 

\\U(J&)-U 5 \\ 2 +aD#(J&,A*) < 5 2 -a(i\A & a -A ] ) = 8 2 -a{{u]\U'{A){A s ol -A))+{£, A s a -A^)). 
By Proposition 8, it follows that 

( W t, U'(A)(A s a - A*)) < (1 + 7 )lMll \\U(A s a ) - U(A) || < (1 + 7)11^11(5+ ||W(^)-W'||). 
Therefore, 

\\U{A 5 a )-U 5 f + aD^{Ai,A ] ) <5 2 + a{l + 1 )\\u^\{5+\\U{Ai)-U s \\) + a\\S\\-\\A 5 a -A ] \\. 

Since ||(f|| is arbitrarily small, it follows that, (£ — ||^||)/C > 0. Moreover, since is 
g-coercive with constant £. 

For the case q = 1, the above inequalities imply that, 

(\\U(A s a ) -U s \\ - q(1 + 7)lMH/2) 2 + a(l - l/CKID^t^^^^ + aCl + ^II^H) 2 

Hence, the assertions follow. 

For the case q > 1, we denote /?i = ||^||/C and we have that, 

q q 

Thus, assuming that (3i = 0(5 1 ' q ), we have that, 

(\\U{A 5 a ) -U s \\ - a i±I|| w t|| > ) + a 1^1 D(i ^AlA ] ) <(6 + a(l + 7)lMll) 2 + a—, 

and the assertions follow. □ 
The rates obtained in Theorem 5 can be seen as a measure of reliability of solutions, 
since it quantifies how better solutions become when the noise level decreases. 
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4 Morozov's Principle 



We establish now a relaxed version of Morozov's discrepancy principle for the specific prob- 
lem under consideration (see [17]). This is one of the most reliable ways of finding the 
regularization parameter a as a function of the data U 5 and the noise level 5. Intuitively, 
the regularized solution should not fit the data more accurately than the noise level. We 
remark that it does not follow immediately because, the parameter now has to be chosen 
as a function of the noise level 5 and the data U 5 . Thus, it is necessary to prove that such 
functional in fact satisfies the required criteria to achieve convergence and convergence rates 
when 5 — » and IA 5 — > U, where here Li represents the noiseless data. 
From Equation (9), it follows that any A G H satisfying 

\\U{A) -IA & \\ < 5 (13) 

could be an approximate solution for (7). If A 5 a is a minimizer of (10), then Morozov's 
discrepancy principle says that the regularization parameter a should be chosen from the 
condition 

\\U(A 5 a )-U s \\=5 (14) 

when it is possible. In other words, the regularized solution should not satisfy the data more 
accurately than up to the noise level. 

In what follows we consider a relaxed condition on Morozov's discrepancy principle, since 
the identity above is too restrictive. See [2]. 

During the following analysis, we also require that the functional /_4 satisfies: 

f Ao (A)=0 A = A . (15) 

Fixing the noise level S and the data U s , we define some auxiliary functionals and sets: 

Definition 5. [2] Define the functionals 

L:AeO. i— >• L(A) = \\U(A)-U 8 \\ GM+U{+oo}, (16) 
H-.AeO. i— ► H{A) = f Ao (A) gR+u{+oo}, (17) 

I:aeR+ _> /(«)=^ a (4)GK + U{+ TO }. (18) 

We also define: 

1. The set of all minimizers of the functional (10) for each a G (0, oo) 

M a := {A s a G : L{a s a ) < L(A), VA e H l (0, S, H 1+£ (D))} . 

Note that we have extended L(A) to be equal to \\U(A) — U s \\ when A G H and be 
equal to +oo otherwise. 

2. The set C of all f ^-minimizing solutions of ( 8) (for this specific problem it is a unitary 
set). 

Definition 6. For 1 < n < r 2 we choose a = a(5,U 6 ) > such that 

Ti5<\\U(A 8 a )-U s \\<T 2 5 (19) 

holds for some A a in M a . 
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We start it by a classical Lemma: 

Lemma 3. [21, Lemma 2.6.1] The functional H(-) is non-increasing and the Junctionals 
!/(•) and I(-) are non- decreasing with respect to a G (0, oo). In other words, for < a < (3 
we have 




Lemma 4. The functional I : (0, oo) — > [0, oo] is continuous. The sets 

M := I a > 



and 

N := la> 



inf L(A d a )< sup L(A d c 



inf H(A s a )< sup H(A S 



are at most countable and coincide. Moreover, the maps L(-) and H(-) are continuous in 
(0,oo)\M. 

Proof: Let {a n } n£ N C (0, oo) be a convergent sequence converging to a*. Thus, we can 
choose a sequence {A n } n£ fq where, for each iigN, A n := A a , a minimizer of the functional 
defined in Equation (10) with a replaced by a n . By the coerciveness of /^ , {^4 n } ne pj is 
bounded and thus, it has a weakly convergent subsequence denoted by {^4fc}fc e N with limit 
A*. 

As U(-) is weakly continuous, the norm of L 2 (0, S, Wl ,2 (D)) is lower semi-continuous and 
fao(-) is weakly lower semi-continuous, we have 

J%A A *) ^ liminf^ iQfc (A) < limsup^, Qfc (A.) 

k— >oo 

which shows that A* G M a *. Thus, by the monotonicity of /(■) we have that /(■) is 
continuous at a*. 

If a G M there are A,Ae M a such that L{A) < L(A), as 1(A) = 1(A), we have 
L(A) ± aH(A) < G(A) ± aH(A) 1(A) - aH(A) < 1(A) - aH(A) <£> H(A) > H(A) 

and M C N. The other inclusion is analogous. The countability of M follows by the fact 
that for each a G M we can associate the interval T a = (mf A 5 aGMa L(A s a ), sup^ eMa L(A s a )). 
By the monotonicity of L we have that for each a, (3 G M, T a R Xp = 0. Therefore, we can 
define an injective map that associates each a G M to an element of X a D Q. 
The continuity of L and H with respect to a out of M follows by the same argument above 
about the continuity of /, with {a n }„ e N and a* in (0, oo)\M. □ 

The following lemma follows from the fact that the sets M a are weakly closed. 
Lemma 5. For each a > 0, there exist A\,A2 G M„ such that 

L(Ax) = inf L(A) and L(A 2 ) = sup L(A) 
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For the remaining part of this Section, we assume that /^ (^4o) = 0. 

Proposition 9. Let 1 < t\ < r 2 be fixed. Suppose that \\U(Ao) — U 5 \\ > t 2 5. Then we can 
find a, a > 0, such that 

L(Ai) < nS < t 2 S < L(A 2 ), 
where we denote Ai : = A s a and A 2 '■= A^. 

Proof: First, let the sequence {a n }„ e N converge to 0. Then, we can find a sequence {A n }neN 
with A n G M an for each n <E N. Now, let A^ be an /^-minimizing solution of (8). Hence, 
we have 

L(A n f < I(a n ) < F% ha M ] ) < 52 + «n/4o(^)- 

Thus, for n sufficiently large, we have L(A n ) 2 < (ti5) 2 , as Oinfj^ia^) — an d f° r this n we 
can set a := a n . 

Now, we assume that a n — > 00 and we take A n as before, the 

H(An) < —I{a n ) < — ?%, an (Ao) = —\\U(A )-U 5 \\ ^ as n -> 00. 

Thus, lim /^ (^4 n ) = 0, which implies that {A n } ne ^ converges weakly to Aq. Then, by the 
weak continuity of U(-) and the lower semi-continuity of the norm, we have 

\\U(Aq) -U 5 \\ < liminf \\U{A n ) -U 5 \\, 

n— >oo 

which shows the existence of a such that L(A~) > t 2 5. □ 
Now we require that, there is no a > with Ai,A 2 G M a such that 

\\U{Ax) - U & \\ < n8 < r 2 5 < \\U(A 2 ) -U s \\. 

If such a exists, A\ would be a sufficiently good approximation for the solution of (8) given 
that ||W(*4.t) — U s \\ < 5. Thus, we can state the following theorem: 

Theorem 6. [2, Theorem 3.10] Under the above assumptions and the assumptions of Propo- 
sition 9, we have the existence of an a := a(5) > and A s a G M a such that 

nS< \\U(A s a )-U s \\ <t 2 5. (21) 

We now present the first main result of this section. It states that, if we use (19) to choose 
a = a(5,U s ), then the regularized solutions A s a converges weakly to the /^-minimizing 
solution of (7). 

Theorem 7. Let 5 > and U 5 satisfy the hypotheses above. Then, the regularizing param- 
eter a = a(5,U s ) obtained through Morozov's discrepancy principle (19) satisfies (11), i.e., 
the limits below 

5 2 

lim a( 8,14 s ) = and lim — — ^ = 0. 

hold. 

Proof: Take a sequence {<5 n }„ S N in (0, +00) converging to zero, fix the noiseless data U and 
denote ot n := a(5 n ,U Sn ) the regularizing parameter chosen from Equation (19). Let A n := 
A 5 a n be a minimizer of (10) with 8 n , a n and U Sn . This defines the sequence {A n } n ^n- As 
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/Uo(') i s coercive, such sequence has a weakly convergent subsequence, denoted by {Ak}km 
with limit A. 

By the weak lower semi-continuity of \\U(-) — U\\ and f Ao , we have the following estimates 
\\U{A)-U\\ < liminf \\U{A k )-U\\ < liminf \\U(A k )-U 5k \\+5 k < liminf(T 2 +l)5 fe = 0, (22) 

ft— >oo ft— >oo ft— >oo 

i.e., A is a solution for the inverse problem. From Equation (19), we have that 

lim \\U{A k ) - U Sk \\ < lim r 2 5 k = (23) 

ft— >oo ft— >oo 

Therefore, for every A^ G C we have the estimate 

r^ 2 k + a k f Ao (Ak) < \\U{A k )-U 5 "\\ 2 + a k f Ao (Ak) < h + a k f Ao (A f ), (24) 
which shows that limsup^^ iU (-4ft) < f Ao (A^). Then, it follows that 

f Ao (A) < lim inf f Ao (A k ) < lim sup f Ao (A k ) < f Ao {A ] ) (25) 

ft— >oo ft— >oo 

for every A^ G C As A^ is an /^-minimizing solution of the inverse problem, we get that 
A G C and f Ao {A) = /^ (^4 t ). This implies that f Ao (A k ) f Ao {A). 

We now prove the first claim in the theorem. Assume that there exists an a > and 
a subsequence {a k } k£ ^ such that a k > a for every k G N. As above select a sequence 
{*4 ft }fteN of minimizers of (10) with 6 k , a k and U Sk . Define further the sequence {*4fc} fte N of 
minimizers of (10) with 5 k , a and U Sk . 

As L in non-decreasing with respect to a and from Equation (19), it follows that 

\\U{A k ) - U &k || < \\U(A k ) - U 5k || < r 2 5 k -> (26) 

On the other hand, 

limsupa/ A) (A) < \imsup(\\U(A k )-U Sk \\ 2 + af Ao (A k )) 

ft— >OD ft— >OQ (07\ 

< ]im S up(\\U(A^)-U Sk \\ 2 + af A ^))=af Ao (A ] ), 1 ' 

ft— >OD 

for every A^ G C, as W(^) = W, the noiseless data. By the coerciveness of f Ao , it follows 
that it has a convergent subsequence, denoted by {^tftj/teN, with limit A G H. Thus, by 
(25), (26), the weak lower semi-continuity of \\U(-) — U\\ and f Ao , it follows that 

\\U{A)-U\\ < liminf || U(A k ) -U\\ < liminf (\\U(A k ) -U Sk \\ + 5 k ) =0 (28) 

ft— >oo ft— >oo 

f Ao (A) < liminf / >4o ( A k ) < lim sup (A k ) < f Ao {A ] ), (29) 



ft— >oo 



ft— >oo 



for every ^ G C. As above, /U (v4.) = f Ao (A^) and thus /^ (w4. ft ) — > f Ao (A^). From the 
above estimates, we have 

\\U(A)-U\\ 2 + af Ao (A) < liminf (\\U(A k ) - U Sk \\ 2 + af Ao (A k )) 

ft— >oo 



< liminf (\\U(A) -U 5k \\ 2 + af Ao {A)) 

ft— >oo 



\\U(A)-U Sk \\ 2 + af Ao (A), 
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for every A G 0, i.e., A is a minimizer for (10) with a and the noiseless data U. Thus, by 
the convexity of f Ao we have, for every t £ [0, 1), 

f Ao ((l - t)A + tA ) < (1 - t)^ (l) + tiU (A)) = (1 - t)f Ao (A). 
On the other hand, from the above estimates it follows that 

af Ao (A) = \\U{A)-Uf + af Ao (A) 

< \\U{{1 -t)A + tA ) -U\\ 2 + af Ao ((l-t)A + tA ) 

< \\U((1 - t)A + tA ) - Uf + 57(1 - t)f Ao (A) 

This implies that atf Ao (A) < \\U((1 - t)A + tA ) - Uf. As U = U(A), by Proposition 5 
with "H = Ao — A, we have 

af Ao (A) < lim -\\U((l-t)A + tA ) -U\\ 2 = 
t->o+ t 

Therefore, f Ao (A) = which, by hypothesis, can only hold if A = Ao- However, by 
hypothesis, for every 5 > 0, Ao is chosen such that ||W(*4.o) ~U S \\ > T\5 > 5. Thus, 

\\U(A )-U\\ > \\U(Ao) -u s \\ - \\U S -U\\ > (n-l)5>0. 

Therefore, we have achieved a contradiction with the fact that U = U(A) = U(Aq). We 
conclude that a(5, U 5 ) — > when 5 — > 0. 

Now for the second limit, take the subsequence {Ak}km of the sequence of the beginning 
of the proof, with 5k — > 0. We know that this sequence converges weakly to an f Ao - 
minimizing solution of the inverse problem A\ with f Ao (Ak) — > f Ao {A^). Thus, from (19), 
it follows that 

rfSl + a k f Ao (A k ) < \\U(A k ) - A 5 * || 2 + a k f Ao (A k ) 
< \\U(Ai)-A s «\\ 2 + UkfAo(^) 
= S 2 k + a k f Ao (A^)- 

This implies that 

□ 

We state now that under the Morozov's choice for the regularization parameter a we 
achieve the same convergence rates of the Theorem 5. This is the second main result of this 
section. 

Theorem 8. Assume that A s a is a minimizer of (10) and a = a(5,U s ) is chosen by Morozov's 
discrepancy principle (19). Then, Theorem 5 holds, i.e., we have the following estimates 

\\U(A s a ) - U(A ] ) || = 0{5) and D e (A s a , A ] ) = 0(5), (30) 
with ^edf Ao (A^)- 

Proof: Assume that A 1 is an /^-minimizing solution of (7) and A s a G M a . The first estimate 
is trivial as 

\\U(A s a ) - U(A) || < \\U{A s a ) -U 5 \\ + \\U & - U{A) || < (r 2 + 1)5. 
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By (19) it follows that 

n5 2 + af Ao (Ai) < \MAl)-U s \\ 2 + af Ao (A 5 a ) 

< \\U(A^-U s r + af Ao (A j ) < 5 2 + af Ao (A j ), 

implying that f Ao (A s a ) < f Ao (A^) as T\ — 1 > 0. Hence, by the definition of Bregman 
distance we have, for every £t e df Ao (A ! ) satisfying the source condition, 

D e (Al^) = fA (A s a )-f Ao (A)-(e,Ai-Ai) 

< 1(^,^-^)1 = \(U>(A^+£,A s a -A<)\ 

< \\^\\\\U'(Ai)(A s a -A)\\ + \\£\\\\A s a -Ai\\ (31) 

< (i + 7)IM||||W(^) -u(A)\\ + ^n^t^L^) 

The last inequality follows by the tangential cone condition and the 1-coerciveness with 
constant £ of f Ao . As $ can be chosen with \\£\\ arbitrarily small, it follows that 

C" 

and then, by (19) 

D e (AlA) < +j)\\J\\\\U(Al) - U(A^)\\ < -r 2 ^— ^j— -(1 + 7 )lM|| • 5. 



□ 



Remark 2. For f Ao q-coercive with q > 1, a reasoning as the one used in Equation (31), 
gives that 

D e {AlA) < ^{AlA^ + m{^i)-U{A^\\ 
^ + -D^(A 5 a ,a) + p 2 \\U(A 5 a )-U(A% 



3 q 

< 

q q 



Assume further that f3\ = 0(Si). As 

\\U{A s a )-U{A)\\ = 0{5) 

it follows that 

\\Ai-A^< l -D !i {AlA^ = 0{5). 
One interesting example is the quadratic functional 

/A)(*4) = M — A) II Hi (CS.-ffi+e (£>))> 

which is 2-coercive with constant 1, since its Bregman distance is 

D2( A - Aq ){*A-,A) = \\A — A\\ H l( 0SH i+ s ( D }y 

Thus, by Theorem 5 it follows that 

0{8) = D 2{A ,_ Ao) (A s a ,A^) = \\A s a - Am Hl(Q ^ H1+HD)y 
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5 Numerical Results 



We first performed tests with synthetic data in order to test accuracy and advantages of the 
methods. Then, we present some reconstructions with real data. 

We note that Problem (2) is solved by a Crank-Nicolson scheme [22]. Since we shall 
use a gradient method to solve numerically the inverse problem, we define the quadratic 
residual J 5 {A) and its gradient V J S (A): 

J* (A) := \\U{A) -U s \\h {0 ,s,mD)) = f \\FM*)) - A*)\\h<p)ds, (32) 



2(U(A) -u s ,w(A)n) LHOtS>LHD)) 

(33) 

{[v(iLyy — iLy)h(t)](s,a(s))}(T,y)drdyds, 



'0 JD 

where, for each s G [0, S], v is the solution of the adjoint equation, 

v T + (av)yy + (av) y + bvy = u(t, a) — u s (s) (34) 

with homogeneous boundary condition. Furthermore, {V : s v(s)} G L 2 (0, S, Wl' 2 {D)). 
We solve Problem (34) numerically also by a Crank-Nicolson scheme. 

For simplicity, in the following examples we assume that the index / of H l (0, S, H 1+e (D)) 
is equal to 1 and considered the standard regularization functional 

fAo(A) = \\A — Ao\\ 2 H i( 0iSH i+e( D }y 



5.1 Examples with Synthetic Data 

Consider the following local volatility surface: 



a(s, u, x) 



0.4(1 - 0.4e- a5 ( u - s )) cos(1.25vr \og(x/s)), if (u, x) G (0, 1.25] x [-0.4, 0.4], 
0.4, otherwise. 

(35) 

The data is generated in a fine mesh and then it is collected in a coarser grid with an 
additive noise. We assume that such noise term is determined by a zero-mean Gaussian 
pseudo-random variable with variance equals to the square of the maximum value reached 
by the noiseless data. We proceed in this way, in order to avoid a so-called inverse crime 
[20]. 

In order to perform the numerical tests, we assume that, r = 0.03, the log-moneyness 
domain is [—5,5] and the time to maturity domain is [0,1]. Moreover, the time and log- 
moneyness step sizes are At = 0.004 and Ay = 0.1, respectively. We also assume that 
s G [29.5,32.5] with different step sizes, As = 0.25,0.1,0.01. 

In what follows, we refer to standard Tikhonov as the case when we consider a single 
price surface. We use the terminology "online" Tikhonov whenever we use more than one 
single price surface. 

Figure 1 illustrates that, when the noise level decreases by increasing the accuracy of 
data, the resulting reconstructions become more similar to the original local volatility sur- 
face, as stated by Theorems 4, 5, 7 and 8. 

In Figure 2, we show that the online Tikhonov presents better solutions than the stan- 
dard one, as we increase the number of price surfaces. Observe that, as it was mentioned 
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above, the parameters are considered in H l (0, S, H l+e (D)) with 1 = 1. Here, the regulariza- 
tion parameter was obtained by the Morozov's discrepancy principle and the regularization 
functional is the standard one. 

Figure 3 shows how the L 2 (D) distance between the reconstructions and the original local 
variance behave as a function of the number of price surface: it is constant for standard 
Tikhonov and non-increasing for online Tikhonov. 




x = log(WSy Time to Maturity ' x=log(Kffi | Time to Maturity x = log(WSy 



Figure 1: Original local volatility surface (image 1) and reconstructions with noise level 
5 = 0.035 (image 2) and 5 = 0.01 (image 3), respectively. As the noise level decreases, the 
reconstructions become more accurate. 




Figure 2: Comparison between standard and online Tikhonov. As the number of price 
surfaces increases, the reconstructions become more accurate. 



5.2 Market Data 

We present now some reconstructions of local volatility by online Tikhonov regularization 
from real data. We make the same assumptions of the previous section about the parameter 
space. We consider seven price surfaces in each experiment. The data corresponds to vanilla 
option prices on futures of WTI oil and Henry Hub natural gas [13]. 

Note that, in order to use the framework developed in the previous sections, we assumed 
that, local volatility is indexed by the unobservable spot price, instead of the future price. 
For more details on such examples, see Chapters 4 and 5 of [1]. 
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Figure 3: L 2 distance between original local variance and its reconstructions, as a function 
of the number of price surfaces, it is constant for standard Tikhonov and non-increasing for 
on line Tikhonov. 

Figures 4 and 5 present such reconstructions for WTI and HH, respectively. We collected 
the data prices for Henry Hub natural gas and WTI oil between 2011/11/16 and 2011/11/25, 
i.e., seven consecutive commercial days. 



WTI Local Vol., case 1 



WTI Local Vol., case 1 
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Figure 4: Local Volatility reconstruction from European vanilla options on futures of WTI 
oil. We used online Tikhonov regularization with the standard functional. 



6 Conclusions 

In this paper we have used convex regularization tools to solve the inverse problem associated 
to Dupire's local volatility model when there is a steady flow of data. We first established 
results concerning existence, stability and convergence of regularized solutions, as conse- 
quences of convex regularization tools and the regularity of the forward operator. We also 
proved some convergence rates. Furthermore, we established a Morozov's discrepancy prin- 
ciple under a general framework, following [2]. Such analysis allowed us to implement the 
algorithms and perform numerical tests. 
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0.55r ; 




Figure 5: Local Volatility reconstruction from European vanilla options on futures of Henry 
Hub natural gas. We used online Tikhonov regularization with the standard functional. 

The main contribution, vis a vis previous works, and in particular [6], is that we extended 
the convex regularization techniques to incorporate the information and data stream that is 
constantly supplied by the market. Furthermore, we have proved a Morozov's discrepancy 
principle that is suitable to this context. 

A natural extension of the current work is the application of these techniques to the 
context of future markets, where the underlying asset is the future price of some financial 
instrument or commodity. In such markets, vanilla options represent a key instrument 
in hedging strategies of companies and in general they are far more liquid than in equity 
markets. The warning here is that, in general, we do not have an entire price surface. 
Actually in this case, we have only an option price curve for each future's maturity. Thus, 
in order to apply the techniques above to this context, it is necessary to assembly all option 
prices for futures on the same instrument (financial or commodity) in a unique surface in 
an appropriate way. This was discussed in [1, Chapter 4] and will be published elsewhere. 
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A Proofs, Technical Results and Definitions 

In this appendix we collect technical results and definitions that were used in the remaining 
parts of the article. We also present the proofs of some results of from Section 2. 

A.l Bregman Distance and g-Coerciveness 

Definition 7. [19, Definition 3.15] Let X denote a Banach space and f : D(f) C X — > 
RUoo be a convex functional with sub-differential df(x) in x G D(f). The Bregman distance 
(or divergence) of f at x G D(f) and £ G df(x) C X* is defined by 

Dt.(x,x) = f(x) - f(x) - (£,x - x), (36) 
for every with (•, ••) the dual product of X* and X . Moreover, the set 

V B {f) = {x G D(f) : df(x)^®} 
is called the Bregman domain of f . 

We stress that the Bregman domain £>b(/) is dense in D(f) and the interior of D(f) is a 
subset of T>B(f). The map x D^(x,x) is convex, non-negative and satisfies D^(x,x) = 0. 
In addition, if / is strictly convex, then D^(x,x) = if and only if x = x. For a survey in 
Bregman distances see [3, Chapter I]. 

Definition 8. For 1 < q < oo and x G D(f), the Bregman distance D^(-,x) is said to be 
q-coercive with constant ( > if 

Ds(y,x)>(\\y-x\\ q x 

for every y G D(f). 

A. 2 Equi- Continuity 

Let X and Y be locally convex spaces. Fix the sets Bx C X and M C C(Bx,Y). A 
set M is called equi-continuous on B x if for every x G B x and every zero neighborhood, 
V C Y there is a zero neighborhood U C X such that G{x$) — G(x) G V for all G G M 
and all x G Bx with x — xq G U. Furthermore, M is called uniformly equi-continuous if 
for every zero neighborhood V C Y there exists a zero neighborhood U C X such that 
G(x) - G(x') G V for all G G M and all x, x' G B x with x - x' G U. 
From [14] we have the technical result: 
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Proposition 10. Let F : [0,T] x Bx — > Y be a function, and Bx, X and Y be as above. 
IfMx := {F(t,-) : t G [0,T]} C C(B X ,Y), M 2 := {F(-,x) : x G B x } C C([0,T],F) and M a 
(respectively M 2 ) is equi- continuous, then F is continuous. Reciprocally, if F is continuous, 
then Mi is equi- continuous and if additionally B x is compact, then M 2 is equi- continuous, 
too. 

A. 3 Proof of Results from Section 2 

Proof of Theorem 1: Well Posedness: Take an arbitrary A G 0, by the continuity of 
A (see Proposition 1) and F, it follows that t i— > F(s,a(s)) is continuous and then weakly 
measurable. Therefore, s i-> \\F(s,a(s))\\ w i,2 {D) is bounded, then U (A) G L 2 (0, S, W 2 (D)), 
which asserts the well-posedness of £/(•). 

Continuity: As F : [0, S] x Q — >■ iy 2 1,2 (.D) is continuous, it follows by Proposition 1 that 
the set {F(s, •) | s G [0, S] } C C(Q, W / 2 1 ' 2 (D)) is uniformly equi-continuous, i.e., given e > 0, 
there is a 5 > such that, for all a, a G Q satisfying \\a — a\\ < 5, we have that 

sup \\F(s, a) — F(s, 5)|| < e. 

sG[0,S] 

Thus, given e > and ^4, ^4 G £3 such that sup s6 r 0iS i \\a(s) — a(s)\\H^(D) < ^ 5 then, by the 
uniform equi-continuity of {F(s, •), s G [0, S]}, it follows that 

\\U{A) - U{A)\\l, {0Aw ^ m = I ~ F(s,~a(s))\\ 2 w ^ {D) ds <e 2 -S, 

J 

which asserts the continuity of U(-). 

Compactness: It is sufficient to prove that, given an e > and a sequence {^4 n }neN C 13 
converging weakly to A, it follows that there exist an uq and a weak zero neighborhood U 
of H l (0, S, H 1+£ (D)) such that for n > n , An - A G U and 

\\U{An) -U{A)\\ L2{0AW i,2 (D)) < e. 

Following the same arguments of the proof of Lemma 1, we can find a set of functionals 
C n ,m G H l (0, S, H 1+e (D))*, defining such zero neighborhood U. We first note that, since 
F is weak continuous, it follows that, given an e > 0, there are a.\, ...,on G H 1+e (D) and 
5 > 0, such that 

sup \\F(s,a) - F(s, a)\\ < e/S 

se[o,s] 

for all a, a G B with 

max{| (a — a, a n )H 1 + e (D) | |n = 1, iV} < & (37) 

By Proposition 1 it follows that (A,a, n )H 1 + e (D) H l ([0,S}) with its norm bounded by 
ll-4||«ll a nll-fl" 1+e (,D)- Then, there is a closed and bounded ball 

A C H l ([0, S]) containing {A, ot n ) h^+^d), for all n — 1, N, 

and A G B. For n = 1, N and the same 5 > of (37), there are f n ,M{n) G ^([0, S 1 ]) 

and £ n > such that, ||/||c([o,s]) < ^ fo r every / G A satisfying 

max{|(/,a n )tfi+e(D) | m = 1, M(n)} < 
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Now, define C n>m := a n <8> f n ,m, with n = 1, N and m = 1, M{n). It is an element of 
H l {0, S, H 1+£ (D))*, where, for each A G H l {0, S, H 1+£ (D)), we have that 

(A, C n ,m)j = ((A, 0! n )iyi+e(£>), fn,m)H l ([0,S]) = + I ' a n) (D)fn,m{k) ■ 

These functionals define a weak zero neighborhood U := r\^ =1 U n with 

E/ n := {A G H l (0, S, H 1+e (D)) : \(A,C n , m )i\ < £„, m = 1, M(n)}. 

Therefore, if {^4.fc}fceN C B converges weakly to A G B, then for a sufficient large fe, Au—A G 
U and by the definition of U, we have that for each n = 1, N, 

> \{A — A,C n>m )l\ = \((A - A, a n )fll-M(D), /n,m)i^([Q,S])| 

for all m = 1, ...,M(n). By the choice of the / n)m G H l ([0, S]), it follows that 

|| (A - A,a n ) H i+s iD) \\ H i( [0iS]) < 5 for all n = 1, ...,iV, 
which implies that ||W(y4.fc) — W(^4)|| L2 ( Q 5iy 1,2 (D)) — e ' 

Weak Continuity: The weak continuity follows directly from the proof of compactness, as we 
use the same framework, only changing the compactness of F, by the weakly equi-continuity 
of {F(s, •) : s G [0, S}} on bounded subsets of Q. 

Weak Closedness: It is true since is weakly closed and U(-) is weakly continuous. □ 
Proof of Proposition 5 By Proposition 4, the family of operators {F(s, ■) : s 6 [0,5]} 
is Frechet equi-differentiable. Take A,U G H l (0, S, H 1+£ (D)), such that A,A + H G £2. 
Thus, define the one sided derivative of W(-) at A in the direction % as 

Z/(.A)ft := {s ^ d a F(s,a(s))h(s)}, 
where for each s G [0, S], dropping t to easy the notation, d a F(s, a)h is the solution of 

— V T + a(Vyy - Vy) + bVy = ll(Uyy ~ Uy) 

with homogeneous boundary conditions and u = u(s,a(s)). From Proposition 3 we have 
the estimate 

\\d a F(s,a{s))h{s)\\ w i,2 {D) < C\\h(s)\\ L 2 {D) \\u yy (s,a(s)) - u y (s,a(s))\\ L 2 (D) . 

Where \\u yy (s,a) — Uy(s,a)\\i?(E>) is uniformly bounded in [0,S] x Q. Thus, U'(A)H is well 
defined and 

S 

< C [ \\h(s)\\ L 2 {D) \\u yy (s,a(s)) - Uy(s,a(s))\\ L 2 (D) ds 
< 



c / \\H s )\\L 2 (D)ds — c\\K\\ H u QtStH1+e , D ^ 
Jo 

Therefore, U'(A) can be extended to a bounded linear operator from H l (0, S, H 1+£ (D)) to 
L a (p,S,W?(D)). 
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Now, let A,H,G € H l (0,S,H 1+£ (D)) be such that, 

a, A + n, A + g,A + n + g e q. 

Denote, v := u(s, a(s) + h(s)) — u(s, a(s)). Thus, 

w := (d a u(s, a(s) + h(s))g(s) - d a u(s, a(s))g{s)) 

satisfies 

-W T + a(Wyy - Wy) = ~ yy ~ Vy] ~ k[(d a U(s , <2 + k) g) yy ~ (<9 a W(S,<2 + k) g) y] , 

with homogeneous boundary conditions (dropping the dependence on s). As above, we have 



u'(A + n)g-u'(A)g " 2 



L*(0,S,W^(D)) J _ vv 2 



\w\\ 2 „ A ,2 lr ^ds 



2 IL. /„ ~/„\\ „. /„ ~/„\M|2 



<ci / \\q(s)\\ L2{D) \\v yy (s,a(s)) - v y (s,d(s))\\ L 2 (D) ds 
°s 

+c 2 J o IIM^IIi-wllM^aW + M*))^)!! W^{D) ds 

<C\\H\\ 

H l (0,S,H 1 + e (D)) \\G\\ 

which yields the Lipschitz condition. □ 
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